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Ph! Abstract 

Q ' We examine the possibility of extracting CP-violating terms in the decay 

,^ ■ K^ — > i:^e^e~ by studying the time evolution of a K^ beam. We focus on the inter- 

^ , ference region and search for clear effects. We find that experiments which average 

'k/ \ over the electron and positron momenta can detect CV violation as an oscillation 

;_, ■ in the decay rate. The branching ratio is (3 — 5) x 10~^^ and direct CV violation 

. . . , dominates over a wide range of the parameters. 
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1 Introduction 

One property of the standard model which is still under active consideration is the origin 
of CV violation. Up to now CP-odd contributions have been observed only in the decays 
of K^ mesons [|l|. In the decay one distinguishes two types of CV violation effects: direct 
CV violation occuring in the amplitudes ( described by e ) and CP-asymmetric terms in 
the mass matrix which is called indirect ( described by e ). The value of e is precisely 
known ( |e| = 2.258 x 10~^ ), but there are still uncertainties concerning the value of e . 
The CERN experiment NA31 found [H 



Re(-) = (2.3 ±0.7) x 10^^ (1) 

while the measurement of the FERMILAB experiment E731 is |^ 

i?e(-) = (0.74 ± 0.59) X 10"^ (2) 

which is still consistent with the predictions of the superweak theory. Thus it is still 
interesting to investigate other processes in order to find if direct CV violation is different 
from zero and providing another crucial test of the standard model. Examples for such 
processes are B meson and rare K'^ meson decays, which are actively investigated both 
theoretically and experimentally. A promising decay channel is K^ —>■ 7T^e~^e~ , where CV 
violation may be relatively large. The specific decay K^ —>■ 7i^e~^e~ has been studied 
extensively and its status was recently reviewed ^ Q. 

We study the time development of this decay channel starting with a pure K^ beam 
and pose the question if one could identify a CP-violating signal in the interference region. 
In particular, we are interested in a signature of direct CV violation which we present in 
this article. We will show that an experiment which studies the time development of K'^ 
decays and averages over the momenta of electron and positron is sensitive to CP- violating 
terms. The new effect appears in the interference region of the Ks and Kl component 
and manifests itself as a time oscillation. The effect follows from general symmetry con- 
siderations as is explained in the next section. In addition, we present an estimate for the 
magnitude of the effect. This experiment is especially suited for laboratories with intense 
K beams like Brookhaven ||^. 

Our paper is organized as follows. In section 2 we derive the formula describing the 
time evolution of a pure K^ state and classify the different contributions. In section 
3 we discuss the calculations available for the amplitudes and their dependence on the 
parameters. Furthermore we give the range of parameters, which is used later on in the 
numerical analysis. In section 4 we present the numerical results for the time evolution 
of the K^ state with special emphasis on CV violation in the interference region. Finally, 
the reder who is interested on the experimental possibility can study section 2 and the 
conclusions in section 4. 



2 Classification of the various amplitudes 

The time evolution of a pure K^ state is given in terms of the time development of the 
physical states Kl and Ks as follows 

\K\t) >= -^ \e-'^^\ \K2 > +e\Ki >) + e-'^''{ \Ki > +e\K2 >)]. (3) 

The decay proceeds through two intermediate states, K^ -^ ir^'-f and K'^ -^ n'^'-f'-f, with 
the single or two photons converting into electron positron pairs. The decaying kaon has 
spin 0, and angular momentum conservation demands the intermediate state n^'-f to be in 
a p-wave. It follows, then, that the CV eigenvalue of 7r°7 is (— l)7r(+l)7(— l)'""*^ = +1. 
Thus Kl can decay through the n^'-f channel in terms of CP-conserving parts, whereas K2 
decays through the CP- violating parts of the Hamiltonian. We denote these amplitudes as 

Ai = < TT^e^e^l H^ \Ki >; CP-conserving ( It gives (4) 

indirect CV violation through e), 
B = <7r°e+e"| H^ |^2 >; CP-violating ( direct ) . (5) 

These would be the only two amplitudes if there were no higher order terms. In fact 
the decay of a i^° can also proceed through the intermediate state 7r°77, which is higher 
order in the electromagnetic coupling. This contribution to the decay width is not a priori 
negligible, because, as we will show, it has to be compared with CP-violating terms which 
are suppressed. To be specific, the intermediate state of a pion and two photons has many 
partial waves so that both CV = +1 and CV = — 1 states are allowed. Thus the decay 

K2 -^ 7r°77 -^ n°e+e- (6) 

is CP-conserving with the final state odd under the CV transformation. In fact the decay 
K2 -^ Tf^'-ff has already been observed. We define the relevant amplitude as 

A2 = < TT^e^e"! H-y-y \K2 >; CP-conserving. (7) 

The decay of a pure K^ beam has the general form 

<7r°e+e-|7^|ir°> (t) = -^ [e-'^^' u{k_)i)K[{B+ + eA^ + A+) + B--i^v{k+) (8) 



V2 



e~^^st -( 



u{k,)i>K[At + e(fi+ + A+) + eB-^^]v{k+)] 



where the -|-, — indicate that the spinors are written out explicitly, with A^*", A2 ^ B^ being 
vector amplitudes and B^ being the axial-vector. For explicit definitions see equations 
(12), (17) and (24). 

The term e{B + A2) is small in comparison to Ai on several resons: 

i) B is small, being CP-violating, 



ii) A2 is small, being higher order in electromagnetism 
iii) these two small terms are multiplied by the small parameter e. 

Neglecting e{B + A2), the Ks decays are CP-conserving. The K^ decays contain A2, which 
is CP-conserving, and the CP- violating amplitudes B and eAi. The amplitude B represents 
direct CV violation, whereas the violation in eAi arises through the mass matrix. Since 
both terms are very likely suppressed, it becomes necessary to consider the A2 term, as 
mentioned above. 

Next, we compute the time evolution of the decays. 



dr 



rr,+r< 



dsdA ^ ' 512 TT^mi 

+ e"^^^ * 2 Re [e"^^™^'* (S+ + A+ + eA+) A+*] } (9) 



X 2 [^(^'"^^'"^')-^' 



where 



S = {Pk- PTvf 

X{s, mj^, ml) = s^ + m% + m^ — 2s m\ — 2s rn^ — 2 m^rn^ and 

This expression shows three time intervals: decays for Ks--, Ki mesons and an interfer- 
ence region. The first two show the typical exponential behavior for the decays, and the 
interference has an oscillatory term as well. We point out an important property of the 
interference term. The A2 amplitude is odd under the CV transformation and thus anti- 
symmetric under the exchange of the electron and positron energies or momenta. This is 
born out by explicit calculation, with equation (^) being linear in A. The Af amplitude 
is even under exchange of the electron and positron momenta. Therefore the term A^A^* 
drops out in an experiment which symmetrizes over the electrons and positrons. The re- 
maining interference terms in equation @ are CP-odd. Thus the presence of an oscillation 
in the interference region is a clear indication of CV violation. In the remaining article we 
estimate each of the amplitudes, we calculate the magnitude of the effect and demonstrate 
it with formulas and several figures. 

3 Estimates for the Amplitudes 

3.1 The direct C7^- violating amplitude B 

Among the amplitudes, B is the best known in the standard model. As discussed by 
several authors, the -B-amplitude is calculated according to Fig. [I| by means of an effective 
Hamiltonian for AS = 1 semileptonic transitions derived by using the operator-product- 
expansion [7—10]. One starts at a high energy scale, where the interaction is point-like, and 



scales down to low energies by means of the renormalization group equation. The procedure 
also includes electromagnetic and strong effects contained in the Wilson coefficients. Since 
the 5-amplitude is CP- violating, it involves the imaginary parts of the Wilson coefficients, 



and the dominant terms are [^ 

ImiCy) O7, and Im{Cs) Og 
with 



O7 



(sL^f^di) {e'j^ e), and Og = (sz, 7^, cii) (e 7^ 756). 



(10) 



(11) 



These Wilson coefficients receive their main contribution from energy scales between rrit 
and rric where perturbative QCD is more reliable. The reduced matrix elements of the 
operators involve quark currents between hadronic states and can be related to Ki^ decays 



through an isospin rotation W^. Neglecting the mass of the electron, the final form of the 
amplitude is 



B 



u{k_)^ 



K 



B' 



B 75 v{k 



(12) 



with 



B' 



.Cjr p 



V^jV;,a2-/2U{3)Im(C,), B' 



Cjrp 

'71 



KjV;>2v^/+(s)/m(Cs). 



Following the Kobayashi Maskawa parametrization for the quark mixing matrix we use for 
the coefficients the values [|10[ 



Im{C7 



1 



VudV. 



- lm{VtdV*) 0.74 and Im{C8 



VudV. 



- lmm,V,:) (-0.70) (13) 



with rrit = 170 GeV. The main uncertainty comes from the factor 

Im{\t) = Im{V^^Vfl) = —siS2S3C2sin6 



(14) 



for which we will allow the range (1.0 — 2.0) x 10 ^. The branching ratio from the B- 
amplitude alone is given by the formula 



dr 

dsdA 



5127r3m|^ 



\B 



+ |2 



\B 



-|2 



A(s,m|-,m^) - A^ 



Varying the parameters, we obtain the range 



BR{Ki 



1^ G a ) direct 



(2.4-9.7) X 10 



-12 



(15) 



(16) 



Later on we will use the -B-amplitude as given above with the corresponding ranges of the 
parameters in order to study the development of a pure K^ beam. 



3.2 The CP-Conserving Amplitude Ai 

The amplitude Ai has the same diagrams ( see Fig. |l| ) as the amphtude B. But the 
approach of using the same effective Hamiltonian as in section ^?l| involves the real parts of 
the Wilson coefficients with large contributions from regions far below rric, where perturba- 
tive QCD is not reliable. For this reason one does not use the QCD effective Hamiltonian, 
but resorts to other low energy methods like chiral perturbation theory. We define Ai 
through the equation 



Ai = < 77° e'^e'in^lKi > = u{k_) AX^k v{k+), 



(17) 



and Af is given by 



^t = ^ K^K; ^8 ^ 2 



W4 



lln^ 
6 m|; 
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with the loop-function 0(s) 



4mi 5 1 

H h - 

3s 18 3 



' Am\ 



1 arctan 1/ \ 



Umj^ 



as calculated in [|T^. Most of the factors here have standard definitions except for gs, 
which is the coupling constant of the octet of pseudoscalar mesons and u^, which is a 
dimensionless coupling constant. Both have to be determined experimentally. From the 
decay K ^ Tin it was found g^ = 5.1. uj+ was determined from a x^-analysis of the 
spectrum for the decay K~^ — > vr+e^e^, based on a calculation of the spectrum in xPT 
|l^ , including the same set of parameters. From their data set the BNL E777 group [jl3| 



derived a value of 



w+ 



0.89 



The decay-width from the Ai amplitude reads 



0.24 
0.14 



(19) 
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Ki 



m^ 



A' 



(20) 



This yields a branching ratio for the decay 



BR{Kl 



7r°e^ 



1.71 X 10"'° - 1.14 X 10 



-12 



(21) 



if the decay occurs only through the CP-violating piece of the K^-Kq mass-matrix. 

In the same way as B, the Ai-amplitude with the corresponding range of parameters 
will be used as an input for the time development of a pure K^ beam. 



3.3 The CP-conserving amplitude A^ 



We have already mentioned that we should include the A2 amplitude even though it is of 
O(a^). Since the decay K^ -^ iT^e'^e~ has not yet been observed it was suggested to study 
the intermediate decay K^ — ^ '7r°77, which has recently been observed, with the branching 
ratio 1 14] 

(1.7 ±0.3) X 10"^ 



BRiKi 



n 



%7) 



22 



Starting from these studies, one should couple the two photons to the final electron-positron 
pair. This is useful but not very direct, because some amplitudes contributing to 
Kl — > 77 are suppressed in the K2 -^ 7r°e~'~e~ amplitude, being proportional to nie and 
their contribution being negligible. As will become clear later on, an amplitude which 
gives a significant contribution to the m^^-distribution gives a very small contribution to 
the semileptonic decay. 

The amplitude for the decay Kl — > 7r°77 has been estimated by two different methods. 



One method uses a two-component model developed by Sehgal and collaborators [jT5|, p!6 
The two contributions included are: (a) a diagram with a charged pion loop to which 
the photons are attached and {j3) vector meson intermediate states. The second method 
applies chiral pertubation theory [|12|. The two approaches differ in several respects, but 
for the amplitude which is dominant in our investigation they agree. This comes about as 
follows: the decay K^ — ^ 7r°77 has several amplitudes, but only one of them is significant 
for A2. It is fortunate that estimates of this amplitude give similar results in the two 



methods. We describe the results of the two-component model |]T^ |18|. 
The amplitude for Kl{pk) — > 7!'^{pn)'y{k)'y{k ) has the general structure 



M 



e^^e" 



F ( k'^K - Qfiuk! -k) 



+ G 



Qtiuk- Pk k' -pk + PKfi Pku k' ■ k 
- PKfi ky k' -Pk - k'^ pku k ■ pk ^ 



(23) 



The pion-loop diagrams contribute only through the F amplitude, whose contribution to 
A2 is proportional to rrie and thus small. This follows by considering the general structure 
of the loop integral and the tensor structure of the term that multiplies the F amplitude 
in equation (^31) . The vector meson pole diagram contributes to the amplitude y4^(s. A) 
defined by 

A2 = A{K2 -^ 7r°e+e-)2^ = u{k^) A+ i>K v{k+). (24) 

Conservation of CV demands that ^2^(5, A) is an odd function of A. Defining 
the absorptive part of A2{s, A) is 



/? = / 



1 



4 mi 



Im Ai 



a Geff A 
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eff 
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- P' 



when P 



llnl±^ 
/5 1-/? 



(25) 



We notice that ImA2 is an odd function of A and in addition the hniit /5 — *> 1 is justified 
for the decay K2 — > 7r°e+e~. In this hmit, equation (|25|) is in agreement with formula (28) 
of ||T^. The formula above has only one coupling Geff which is chosen in such a way that 
it reproduces the branching ratio in equation (^21). 



The dispersive part is calculated with the help of a dispersion relation |T8 

ImA+(A) 



Re A^ 
where the lower limit is given by 



TT 



dz 



z — s 



(26) 



^m,zn(A) 



Ami 



-\ -1 



1 - 



{{rriK + Tn^^Y ~ Ami) {{^^k — "^tt)^ ~ AwlI) 



(27) 



and the upper limit is determined by the heaviest particle considered, this being the p 
meson. The differential decay-width reads 



dV 

dsdA 



51277^771^ 



1^2 I 



K 



A(s,m 



2 



m^ 



-A' 



(2J 



Inserting the numerical values we obtain a branching ratio dominated by the vector meson 
coupling constant 



BriKi 



7r°e^ 



l2j 



AM X 10 



-12 



G 



eff^K 



with p 



and G, 



eff^K 



4.15 X 10"^2 

^ displ J- abs 

0.15 X 10"^ 



0.25 X 10 



1.5 



-7 



:i+p) 



(29) 



(30) 



The second method for calculating the decay Kl —>■ 77^77 is chiral perturbation theory. 
The authors include effects of the order p^ ^0\ and p^ [21| in the momentum expansion of 
xPT, as well as vector mesons p2i. This enables them to reproduce the observed decay 



rate and spectrum. In this approach there are again two effective coupling constants which 
have to be fixed experimentally. Here the vector meson coupling constant was chosen in 
such a way so that the measured decay rate is reproduced. 

The calculation of the two-photon-exchange contribution to the decay Kl — > 7T^e~^e~ on 
the basis of the xPT prediction is analogous. The branching ratio achieved in this manner 
is dominated by the vector meson intermediate state and reads 



BriKi 



7r°e+e- 



J27 



1.8 X 10 
4.5 X 10 



-12/ 



-12 



P) 



(31) 



with the same p. In the numerical analysis of section ^ we can choose any of the two 
calculations for K2 -^ 7r°e+e~|27, because the term relevant for our purpose is practically 
the same. In addition, we shall demonstrate that in experiments which average over the 
e"*" and e~ momenta, this A2 amplitude drops out in the interference region. 



8 



4 Numerical results 

With the amphtudes developed and the general equation (H) we calculate the time evolution 
for the decay K'^ — > 7r°e"'"e~. As already discussed, there are three regions of physical 
interest . 

1) The i^s-region, where the CP-conserving amplitude contributes to the decay width. 
A measurement in this region will determine parameters of the Ks-decaj like u.^- . 

2) The i^L-decay region, which has been studied in several articles interested in CV 
phenomena ( see the review [Q |^ )• The relevant amplitude in this case is 15 + ^42 + 
e^ip, which gives several terms. The interference term Re{A2B^ + €^1^2) is odd in 
A and one could define an asymmetry in A in order to extract this term. The sum 
of the absolute values squared is even in A, and we will need precise measurements 
of the amplitudes in order to observe an excess of events. 

3) More interesting is the interference region occuring in the time interval (8 — 9) • r^-^. 
This term has an oscillatory behavior. The term A2A1 is a linear function in A and 
drops out when we average over the electron positron pair. The remaining terms 
BAi and ejAip are both CP-violating. Thus the appearance of an oscillation in the 
interference region gives evidence for CV violation. 

We have studied this phenomena numerically and show the effect in several figures. 
Fig. (2) and Fig. (3) show the branching ratio as a function of time for two different 
time scales. We note that an oscillation is evident. We use u;+ = 0.89 and three values 
for ImXt = 1.0 X 10~^, 1.5 x 10~^ and 2.0 x 10~^. We plotted the same curves in 
figures (4) and (5) where the CP- violating terms BAf and e|v4ip are set equal to zero. 
We noticed that the curves to be compared are different. In the interference region 
there is a clear oscillation and for very long times the curves which contain the CV 
amplitudes lie above the curve without the CP-violating terms. For the latter region 
experiments studying i^^-decays need a precise measurement of magnitudes in order 
to establish a signal. The branching ratio is in the range (3 — 5) x 10~^^. In contrast to 
this situation, the oscillation in the interference region is unambiguous. A comparison 
of the magnitudes of the contributing terms shows that BAi dominates over e|Aip. 
Only if we choose uj+ at the upper bound (u;+ = 1.13) the two contributions are of 
comparable size. But the direct CP- violating term is still larger by (6 — 113)% for 
Im{Xt) = (1.0-2.0) X 10"^ 

We conclude that an experiment searching for a branching ratio down to 10~^^ and sensitive 
to the time development of the decay can observe CV violation as an oscillation in the 
interference region. The experiment does not require a measurement of the e~^e~ energy 
asymmetry. 
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Figure Captions 

Fig. 1^ Short distance contributions to the decay K — > 7r°e+e~ 

Fig. 2 Time development of the partial branching ratio r(i^'^ -^ 'K^e^e^){t)/T{KL -^ all) 
for the 
time interval (6 — 20) ■ t^s- 

Fig. 3 Same as in fig. 2 for a larger time interval. 

Fig. 4 The same curves as in fig. 2 together with the modified decay rate in which the 
CP-violating terms are set to zero. 

Fig. 5 The same curves as in fig. 3 together with the modified decay rate when the CV- 
violating terms are set to zero. 
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Figure 1: Short distance contributions to the decay K -^ Ti^e^e 
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Figure 2: Time development of the decay r{K^ -^ vr^e+e /r{K^ -^ all) for the time 
interval (6 — 20) r^-g 
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Figure 3: Same as in fig. 2 for a longer time interval. 
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Figure 4: The same curves as in fig. 2 together with the decay when the CP-violating 
terms are set to zero. 
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Figure 5: The same curves as in fig. 3 together with the decay when the CP-violatin^ 
terms are set to zero. 
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